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We study the ground-state properties of the S = 1/2 Heisenberg models on the quasi-one- 
dimensional kagome strip lattices by the exact diagonalization and density matrix renormalization 
group methods. The models with two different strip widths share the same lattice structure in their 
inner part with the spatially anisotropic two-dimensional kagome lattice. When there is no mag- 
netic frustration, the well-known Lieb-Mattis ferrimagnetic state is realized in both models. When 
the strength of magnetic frustration is increased, on the other hand, the Lieb-Mattis-type ferrimag- 
netism is collapsed. We find that there exists a non-Lieb-Mattis ferrimagnetic state between the 
Lieb-Mattis ferrimagnetic state and the nonmagnetic ground state. The local magnetization clearly 
shows an incommensurate modulation with long-distance periodicity in the non-Lieb-Mattis ferri- 
magnetic state. The intermediate non-Lieb-Mattis ferrimagnetic state occurs irrespective of strip 
width, which suggests that the intermediate phase of the two-dimensional kagome lattice is also the 
non-Lieb-Mattis-type ferrimagnetism. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.40.Mg 
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I. INTRODUCTION 

Ferrimagnetism is a fundamental phenomenon in the 
field of magnetism. The simplest case of the conven- 
tional ferrimagnetism is the ground state of the mixed 
spin chain. For example, there is an (s, S)=(l/2, 1) 
mixed spin chain with a nearest-neighbor antiferromag- 
netic (AF) isotropic interaction. In this system, the 
so-called Lieb-Mattis (LM)-type ferrimagnetismir— is re- 
alized in the ground state because two different spins 
are arranged alternately in a line owing to the AF in- 
teraction. Since this type of ferrimagnetism has been 
studied extensively, the magnetic properties and the oc- 
currence mechanism of the LM ferrimagnetism are well 
known. In particular, the ferrimagnetism in the quan- 
tum Heisenberg spin model on the bipartite lattice with- 
out frustration is well understood within the Marshall- 
Lieb-Mattis (MLM) theorem^. In the LM ferrimagnetic 
ground state, the spontaneous magnetization occurs and 
the magnitude is fixed to a simple fraction of the satu- 
rated magnetization. 

In recent years, a new type of ferrimagnetism, which is 
clearly different from LM ferrimagnetism, has been found 
in the ground state of several one-dimensional frustrated 
Heisenberg spin systems 7 - - — . The spontaneous magneti- 
zation in this new type of ferrimagnetism changes grad- 
ually with respect to the strength of frustration. The 
incommensurate modulation with long-distance period- 
icity in local magnetizations is also a characteristic be- 
havior of the new type of ferrimagnetism. Hereafter, we 
call the new type of ferrimagnetism the non-Lieb-Mattis 
(NLM) type. The mechanism of the occurrence of the 
NLM ferrimagnetism has not yet been clarified. 

On the other hand, some candidates of the NLM ferri- 
magnetism among the two-dimensional (2D) systems, for 
example, the mixed-spin J\-Ji Heisenberg model on the 
square lattice^ and the 3 =1/2 Heisenberg model on the 



Union Jack lattice^, were reported. These 2D frustrated 
systems have the intermediate ground state, namely, the 
"canted state" , in which the spontaneous magnetization 
is changed when the inner interaction of the system is 
varied. It has not been, however, determined whether 
the incommensurate modulation with long-distance pe- 
riodicity exists in the local magnetization of the canted 
state owing to the difficulty of treating these 2D frus- 
trated systems numerically and theoretically. Therefore, 
the relationships between the canted states of these 2D 
frustrated systems and the NLM ferrimagnetic state are 
still unclear. 

Under these circumstances, recently, another candi- 
date of the NLM ferrimagnetism among the 2D systems 
was reported in ref . I VII in which the S — 1/2 Heisen- 
berg model on the spatially anisotropic kagome lattice 
depicted in Fig. Ufa) was studied. A region of the 
intermediate-magnetization states is observed between 
the LM ferrimagnetism that is rigorously proved by the 
MLM theorem^ and the nonmagnetic state of the spa- 
tially isotropic kagome-lattice antifcrromagnet in the ab- 
sence of magnetic fielcU£~— . The local magnetization in 
the intermediate state of the kagome lattice was investi- 
gated by the exact diagonalization method, and it was 
reported that the local magnetization greatly depends 
on the position of the sites, although it is difficult to 
determine clearly whether the incommensurate modula- 
tion with long-distance periodicity is present. This result 
leads to the high expectation that the intermediate state 
of the spatially anisotropic kagome lattice is the NLM 
ferrimagnetic state. Additional research is desirable to 
conclude that the intermediate state of this 2D system is 
the NLM ferrimagnetism. 

In this paper, we study the ground-state properties 
of the S — 1/2 Heisenberg models on the quasi-one- 
dimensional (Q1D) kagome strip lattices depicted in 
Figs. Hfb) and Q^c) instead of the 2D lattice depicted 
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FIG. 1: Structures of the lattices: the spatially anisotropic kagome lattice (a) and the quasi-one-dimensional kagome strip 
lattices (b) and (c) with different widths. An S = 1/2 spin is located at each site denoted by a black circle. Antiferromagnetic 
bonds Ji (bold straight line) and J2 (dashed line), and ferromagnetic bond Jf (dotted line). The sublattices in a unit cell of 
lattice (b) are represented by A, A', B, B', C, C, and D. The sublattices in a unit cell of lattice (c) are represented by A, A', 
B, and C. 



in Fig. (Ha). Note that the inner parts of the lattices in 
Figs.QTb) and [He) are common to a part of the 2D lat- 
tice in Fig. QJa) . We also note that the lattice shapes of 
strips in the present study are different from some kagome 
strips (chains) studied in refs. IVIIIVI1 where the nontriv- 
ial properties of kagome antiferromagnets were reported. 

According to the study in ref. IVI1 it was already known 
that the NLM ferrimagnetism is realized in the ground 
state of the kagome strip lattice in Fig. 0Jc). In the 
present study, we show that both the lattice in Fig. [TJc) 
and the lattice in Fig. [ljb) reveal the NLM ferrimag- 
netism in the ground state. Note also that the lattice 
shape at the edge under the open boundary condition 
depicted in Fig. [TJc) is different from that in ref. I VII [see 
Fig. QJb) in ref. I VI] . Thus, one can recognize that the 
results of the strip lattice with small width are irrespec- 
tive of boundary conditions. We also present clearly the 
existence of the incommensurate modulation with long- 
distance periodicity in the local magnetizations of both 
models in Figs.QJb) and[Uc). Our numerical calculations 
suggest that the intermediate state of the 2D lattice in 
Fig. [TJa) is the NLM ferrimagnetism. 

This paper is organized as follows. In §2, we first 
present our numerical calculation methods. In §3, we 
show the ground-state properties of the lattice depicted 
in Fig. [TJc) in finite-size clusters. In §4, we show the 
ground-state properties of the lattice depicted in Fig. 
QJb). Sections 5 and 6 are devoted to discussion and 
summary, respectively. 



renormalization group (DMRG) metho d 32 ' 33 . 

The ED method is used to obtain precise physical 
quantities of finite-size clusters. This method does not 
suffer from the limitation of the cluster shape. It is ap- 
plicable even to systems with frustration, in contrast to 
the quantum Monte Carlo (QMC) method coming across 
the so-called negative-sign problem for systems with frus- 
tration. The disadvantage of the ED method is the lim- 
itation that available system sizes are very small. Thus, 
we should pay careful attention to finite-size effects in 
quantities obtained by this method. 

On the other hand, the DMRG method is very powerful 
when a system is (quasi-)one-dimensional under the open 
boundary condition. The method can treat much larger 
systems than the ED method. Note that the applicability 
of the DMRG method is irrespective of whether or not 
the systems include frustrations. In the present research, 
we use the "finite-system" DMRG method. 



III. KAGOME STRIP LATTICE WITH SMALL 
WIDTH 



II. NUMERICAL METHODS 

We employ two reliable numerical methods: the ex- 
act diagonalization (ED) method and the density matrix 



In this section, we study the magnetic properties in 
the ground state of the 5 = 1/2 Hcisenberg model on the 
kagome strip lattice depicted in Fig. [TJc). The Hamilto- 
nian of this model is given by 
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FIG. 2: (Color) (a) Dependence of the lowest energy on Stot- 
The results of J2/J1 = 0.5 and 0.58 for the system size of 
N = 96 are presented. Arrows indicate the values of the 
spontaneous magnetization M in each J2/J1. The position 
of an arrow is given by the highest S£ ot value among those 
that give the common lowest energy, (b) J2/J1 dependence 
of M/M B obtained from ED calculations for N = 24 (black 
square) under the periodic boundary condition and DMRG 
calculation for N = 48 (red triangle) and 96 (blue inverted 
triangle) under the open boundary condition. Note that for 
N = 96 in (a) and (b), we use MS = 900 and SW = 15 and, 
that for N = 48 in (b), we use MS = 400 and SW = 15. 



where Sj^ is the S — 1/2 spin operator at the £-sublattice 
site in the i-th unit cell. The positions of the four sub- 
lattices in a unit cell are denoted by A, A', B, and C in 
Fig. HJc). Energies are measured in units of Ji ; we fixed 
J 1 = 1 hereafter. In what follows, we examine the region 
of < J2/J1 < 00 in the case of Jf = — 1. Note that 
the number of total spin sites is denoted by N; thus, the 
number of unit cells is N/4. 

We examine the J2/J1 dependence of the ratio M/M s , 
where M and M s are the spontaneous and saturation 
magnetizations, respectively. Let us explain the method 
used to determine M as a function of N and Jij J\. First, 
we calculate the lowest energy E(Jij 'Ji, S* ot , N), where 
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FIG. 3: (Color) Local magnetization {S*e) at each sublattice 
£; A (black square), A' (red triangle), B (blue square), and C 
(purple inverted triangle). Panels (a) and (b) show results for 
J2/J1 = 0.3 and 0.57, respectively. These results are obtained 



from our DMRG calculations for N = 96 (i 



■,24). 



S* ot value is the z-component of the total spin. For ex- 
ample, the energies for various S* ot in the two cases of 
J2/J1 are presented in Fig. EJa); the results are obtained 
by our DMRG calculations of the system of N = 96 with 
the maximum number of retained states (MS) of 600 
and the number of sweeps (SW) of 15. The spontaneous 
magnetization M(J 2 / J\, N) is determined as the highest 
St ot among those at the lowest common energy [see ar- 
rows in Fig. Ufa)]. Our results of the J2/J1 dependence 
of M/M s are depicted in Fig. [U[b). We find the existence 
of the intermediate magnetic phase of < M/M s < 1/2 
between the LM ferrimagnetic phase of M/M a = 1/2 
and the nonmagnetic phase. In order to determine the 
spin state of this intermediate phase, we calculate the 
local magnetization (<5?c), where (A) denotes the expec- 
tation value of the physical quantity A and S^ is the 
z-component of Sj^. Figure [3] depicts our results for a 
system size N — 96 on the lattice depicted in Fig. [TJc) 
under the open boundary condition; Figs. E^a) and|3|b) 
correspond to the case of the LM ferrimagnetic phase and 
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that of the intermediate phase, respectively. 

The results clearly indicate the existence of the incom- 
mensurate modulation with long-distance periodicity in 
the intermediate phase. We also confirm that the period- 
icities of the local magnetizations in the NLM ferrimag- 
netism in the present model depend on the J2/J1 value 
but not on the length of the system, as reported in the 
case of ref. IVII 

It is worth emphasizing the point that the intermediate 
phase commonly exists irrespective of the shape at the 
edge of the strip lattice when one compares the result 
of the present lattice depicted in Fig. [TJc) and that in 
ref. IVII Therefore, we can conclude that the intermediate 
phase exists and that it is NLM ferrimagnetic. 

Although there exists an intermediate NLM ferrimag- 
netic phase in the case of the kagome strip lattice de- 
picted in Fig. HJc), we should note that there is a large 
discrepancy in dimensionality between the kagome strip 
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Here, the positions of seven sublattices are denoted by A, 
A', B, B', C, C, and D in Fig. QJb). Note that the last 
term of eq. [5] is the Zeeman term. The number of spin 
sites is denoted by N. The number of unit cells is N/7; we 
consider N/14 as an integer. We mainly use the DMRG 
method for investigating the magnetic properties in the 
ground state of this Q1D system under the open bound- 
ary condition. We also investigate the properties under 
the periodic boundary condition by the ED method, al- 
though the size treated by this method is only in the case 
of N = 28. Hereafter, we consider J\ = 1 as an energy 
scale and we investigate the region of < J2/J1 < 00 in 
the case of Jf = — 1. 



B. Phase diagram 

First, let us examine the J2/J1 dependence of the ratio 
M/M s in the absence of the external magnetic field h. 
The procedure for determining M is the same as that 
mentioned in §3 [see also Fig. Ufa)]. We present our 
results of the spontaneous magnetization in Fig. HJb). 
We successfully observe the intermediate-magnetization 



lattice depicted in Fig. [ljc) and the 2D kagome lattice 
depicted in Fig. HJa). In the next section, we treat the 
kagome strip lattice depicted in Fig. [ljb) whose width is 
larger than that of the kagome strip lattice depicted in 
Fig. W)- 



IV. KAGOME STRIP LATTICE WITH LARGE 
WIDTH 

A. Hamiltonian 

In this section, we study the ground-state properties 
of the S — 1/2 Heisenberg model on the kagome strip 
lattice depicted in Fig. [IJb). The Hamiltonian of this 
model is given by 
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region irrespective of the boundary conditions. A careful 
observation of Fig. Hfb) enables us to observe the eight 
regions at least in the finite-size system. As a matter 
of convenience, hereafter, we call these regions Ri , R2 ? 
■ •-, R7, and Rg. In the case of N — 112 under the open 
boundary condition, for example, Fig. [5ja) illustrates the 
regions Ri to R 8 : Ri is the region of M/M s — 3/7, R2 is 
the region of 11/28 < M/M s < 3/7, R 3 is the region of 
1/8 < M /M s < 11/28, R 4 is the region of M/M s = 1/8, 
R 5 is the region of < M/M s < 1/8, R6 is the region of 
M/M s = 0, R 7 is the region of < M/M s < 1/7, and R 8 
is the region of M/M s = 1/7. Here, the dashed lines in 
Fig. [3a) indicate the boundaries of these regions. 

It should be noted that the values of M /M s in the 
R4 region and that at the lower edge of the R2 region 
change with increasing N, as shown in Fig. @|b); the for- 
mer value is M/M s = (N — IA)/7N and the latter value 
is M/M s = (3N - 28)/7A. These changes due to the in- 
crease in system size come from the finite-size effect. We 
find that the value of M/M s in the R4 region under the 
open boundary condition increases and approaches the 
value of M/M s = 1/7 when N increases. Furthermore, 
the magnetization value in the R4 region is M/M s = 1/7 
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FIG. 4: (Color) (a) Lowest energy in each subspace divided by 
Stot- The results of the DMRG calculations obtained when 
the system size is N = 56 for J2/J1 = 0.20 and 0.57 are 
presented. Arrows indicate the spontaneous magnetization M 
for a given J2/J1. (b) J2/J1 dependence of M/M B obtained 
from the ED calculations for N =28 (blue triangle) under the 
periodic boundary condition and the DMRG calculations for 
N =56 (red square) and 112 (black pentagon) under the open 
boundary condition. Note that for N = 56 in (a) and (b), 
we use MS = 600 and SW =15, and that for N = 112 in 
(b), we use MS > 900 and SW =15. Here, MS denotes the 
maximum number of retained states and SW the number of 
sweeps used in DMRG calculations. 



in the case of N — 28 under the periodic boundary condi- 
tion. Therefore, it is expected that the value of M/M s in 
the R4 region is 1 /7 in the thermodynamic limit. We also 
confirm that the value of M/M s at the lower edge of the 
R2 region gradually increases and approaches the value 
of M/M s = 3/7 with increasing N. In addition, we can- 
not confirm the R2 region in the case of N — 28 under 
the periodic boundary condition. These circumstances 
indicate a possibility that the R2 region merges with the 
Ri region of M/M s = 3/7 in the thermodynamic limit. 

Next, to determine whether each region survives in 
the thermodynamic limit, we study the size depen- 
dences of the boundaries between the regions under the 
open boundary condition. Figure |5fb) shows the re- 
sults of N =42, 56, 84, and 112 from DMRG calcu- 
lations. Note here that we define R2 as the region of 



(3N - 28)/7N < M/M s < 3/7 and R 4 as the region of 
M/M s = (N - 14) /7N in the finite-size system. One 
can find immediately from Fig. [5jb) that all regions, ex- 
cept the R7 region, survive in the limit N — > 00. To 
determine whether the R7 region survives in the thermo- 
dynamic limit, we investigate the size dependence of the 
width of the R7 region in Fig. [Tj] This plot shows us that 
the width of the R7 region decreases with increasing N. 
It is difficult to determine whether the R7 region survives 
in the thermodynamic limit. The convex downward be- 
havior is observed for large sizes so that the region might 
survive; however, the observed behavior may be one of 
the serious finite-size effects. The issue of establishing the 
presence or absence of the R7 region should be clarified 
in future studies. 
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FIG. 5: (a) Definitions of the Ri-Rs regions in the case of 
N = 112 under the open boundary condition: Ri is the region 
of M/M s = 3/7, R 2 is the region of 11/28 < M/M s < 3/7, 
R 3 is the region of 1/8 < M/M s < 11/28, R 4 is the region of 
M/M s = 1/8, R 5 is the region of < M/M s < 1/8, R 6 is the 
region of M/M s = 0, R 7 is the region of < M/M s < 1/7, 
and Rg is the region of M/M s = 1/7. Dashed lines indi- 
cate the boundaries of these regions, (b) Size dependences 
of boundaries under the open boundary condition. The re- 
sults presented are those of N =42, 56, 84, and 112 from 
DMRG calculations. The curve line named R;-R;+i indicates 
the boundary line between the R; and R;+i regions, where I 
is integer. 
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FIG. 6: Size dependence of the width of the R7 region. It is 
difficult to determine whether the R7 region survives in the 
thermodynamic limit. 



C. Magnetic properties in each region 

In this subsection, we investigate the local magnetiza- 
tion (Sf^) to study the magnetic structures in various 
regions except the R2 and R6 regions. Note that we cal- 
culate (Sf^) within the subspace of the highest S* ot cor- 
responding to the spontaneous magnetization M when 
J2/J1 is given. Considering the fact that the present 
lattice depicted in Fig. QJb) has seven sublattices in the 
system, we will use different colors or symbols for each 
sublattice £ for presenting our results of (Sf^), as de- 
picted in the inset of Fig. [TJa); we use a black square for 
£ =A, a red triangle for £ = A', a blue cross for £ = B, a 
green pentagon for £ = B', a purple inverted triangle for 
£ =C, an aqua diamond for £ = C, and a black circle for 

e=D. 

First, we examine the Ri and Rs regions. We present 
our DMRG results of (Sf^) of the system of N = 112 in 
Figs. [Tfa) andJTJb) for J2/J1 — 0.2 and 1.9, respectively. 
In Fig. [TJa) , we observe the uniform behavior of upward- 
direction spins at the sublattice sites A, A', B, B', and 
D, and downward-direction spins at the sublattice sites 
C and C. In Fig. [TJb), we also observe the uniform 
behavior of upward-direction spins at the sublattice sites 
B, B', C, and C, and downward-direction spins at the 
sublattice sites A, A', and D. Therefore, we conclude that 
the LM ferrimagnetic states are realized in the regions of 
Ri and Rg. 

Our understanding of the origins of these LM ferri- 
magnetic phases is based on the Marshall-Lieb-Mattis 
(MLM) theorem. In the case of J2/J1 = 0, no frustra- 
tion occurs; thus, the spin state depicted in Fig. Ufa) is 
realized. This state shows the LM ferrimagnetism with 
M /M s — 3/7. The Ri region is directly connected to the 
LM ferrimagnetic state of J2/J1 = 0. Therefore, the Ri 
region of M/M s — 3/7 is regarded as the LM ferrimag- 
netic phase. In the limit J2 — > 00, on the other hand, the 
present model becomes equal to a model of an S = 1/2 
diamond chain depicted in Fig. [5Jb). The value of mag- 
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FIG. 7: (Color) Local magnetization (Sfe) at each sublattice 
£. Panels (a) and (b) show results for J2/J1 =0.2 and 1.9, 
respectively. These results are obtained from our DMRG cal- 
culations for TV = 112 [i =1, 2, • ■ •, 16). The inset in the 
panel (a) presents the correspondence relationship between 
each colored symbol and each sublattice £ used in Figs. [7] [9] 
and [10] 



netization takes M/M s = 1/7 in the ground state of this 
diamond chain according to the MLM theorem. There- 
fore, the Rs region of M/M s = 1/7 is regarded as the LM 
ferrimagnetic phase. 

Next, we investigate the R3, R5, and R7 regions. Our 
results obtained from the DMRG calculations of N = 168 
are depicted in Figs. M,&)W C ) for J2/J1 = 0.57, 1.14, 
and 1.69, respectively. We clearly observe incommensu- 
rate modulations with long-distance periodicities in every 
case in Fig. ® In addition, we confirm from Fig. |U[b) 
that the ratio M/M s changes gradually with the varia- 
tion in J2/J1 in the R3, R5, and R7 regions. Since the 
widths of the R3 and R5 regions survive in the thermo- 
dynamic limit as was clarified in the previous subsection, 
we conclude that the R3 and R5 regions are NLM ferri- 
magnetic phases. Although it is unclear whether the R7 
region survives in the thermodynamic limit, this region 
is an NLM ferrimagnetic phase if it survives. 

Finally, in this subsection, we examine the R4 region. 
Our result of (Sf^) for J2/J1 = 1 in the system of 
N = 168 is depicted in Fig. QJJ] We do not detect the in- 
commensurate modulation in this R4 region. In addition, 
we confirm from Fig. 0Jb) that the ratio M/M s in the 
R4 region does not change with the variation in J2/J1 in 
contrast to the cases in the R3 and R5 regions. These 
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FIG. 8: (a) Kagome strip lattice depicted in Fig. \V[b) in the 
limit of J2/J1 = 0. White arrowheads denote the classical 
spin configuration in the LM ferrimagnetic state of M/M s = 
3/7. (b) Kagome strip lattice depicted in Fig. QJb) in the 
limit of J2 — > 00. White circles represent the effective S = 
1/2 spins formed by five S = 1/2 spins at the sublattices 
£ =A, A', B, B', and D. The black bold and black dotted lines 
show the antiferromagnetic and ferromagnetic interactions, 
respectively. 



circumstances suggest that the R4 region is the LM ferri- 
magnetic phase. However, one cannot speculate the spin 
state from the result of (S^) because it is difficult to 
clarify the spin state on the basis of the results of deter- 
mining whether the spins are up or down, as was success- 
fully observed in the Rs region. We further investigate 
the magnetization curve in this region to know whether 
the R4 region is the LM ferrimagnetic phase. The mag- 
netization curves of J2/J1 = 1 for N = 56 and N = 112 
calculated by the DMRG method are presented in Fig. 
nUa) 34 . Figure \TV[b) is obtained by zooming the region 
of h near h = 0. One can observe the existence of the 
magnetization plateaus at the height of the spontaneous 
magnetization for both system sizes. We also confirm 
that the difference in the width of the plateau between 
N = 112 and 56 is very small. These features indicate 
that the spin gap exists in the R4 region in the thermo- 
dynamic limit. If the R4 region is the NLM ferrimagnetic 
phase, no spin gap is present in this region. (It was re- 
ported in ref. IVII that the NLM ferrimagnetism is gapless 
as a response to a uniform magnetic field.) Therefore, 
we conclude that the R4 region is the LM ferrimagnetic 
phase. 

V. DISCUSSION 

We discuss the relationships between the interval < 
J2/J1 < 1 of the kagome strip lattice depicted in Fig. 
QJb) and those of the spatially anisotropic kagome lat- 
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FIG. 9: (Color) Local magnetization (Sf^) at each sublat- 
tice f. The correspondence relationship between each colored 
symbol and each sublattice £ is described in the inset in Fig. 
[7ja). Panels (a), (b), and (c) show results for J2/J1 =0.57, 
1.14, and 1.69, respectively. These results are obtained from 
our DMRG calculations for TV = 168 (i =1, 2, • • •, 24). 



tice studied in ref. I VII while we consider the case of 
another new lattice, with a larger but finite width along 
the direction perpendicular to the bonds of interaction 
Jp in Fig. Gib), namely, the strip width. 

The ratio M/M s of the Ri and R2 regions is commonly 
M/M s — 3/7 in the thermodynamic limit of the lattice 
depicted in Fig. rrfb). The difference of M/M s = 3/7 
from M/M s = 1/3 in the case of the LM ferrimagnetic 
phase of the spatially anisotropic kagome lattice in ref. 
IVIl is attributed to the finiteness of the strip width. Thus, 
the ratio M/M s approaches M/M s = 1/3 when the strip 
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FIG. 10: (Color) Local magnetization (SJt) at each sublattice 
£. This result is obtained from our DMRG calculations for 
J2/J1 = 1 in the system of N = 168 (i =1, 2, • • •, 24). The 
correspondence relationship between each colored symbol and 
each sublattice £ is described in the inset in Fig. [TJa). 



width increases; the Ri and R2 regions of the present 
model depicted in Fig. [lib) correspond to the LM ferri- 
magnetic phase of the spatially anisotropic kagome lat- 
tice studied in ref. IVIl 

The ratio M/M s = 1/7 at J2/J1 = 1 in the present 
model depicted in Fig. \V[b) may be related to the fact 
that the model includes seven sublattices, namely, the 
strip width is finite. At least at J2/J1 — 1, on the other 
hand, it is widely believed that the spontaneous mag- 
netization disappears in the limit of infinite widthi^r— . 
Although the relationship should be clarified in the ex- 
amination of the case of the lattice with an even larger 
strip width, there are two possibilities of the behavior 
near the case of J%] J\ — 1. One is the case when the 
ratio M/M s decreases with increasing strip width and fi- 
nally vanishes, while the LM ferrimagnetic phase, such 
as R4, survives in systems with larger strip widths. In 
this case, the R3 region corresponds to the NLM phase 
of the two-dimensional model on the spatially anisotropic 
kagome lattice. In the other case, the LM ferrimagnetic 
phase, such as R4, becomes narrower with increasing 
strip width, while the ratio M/M s does not vanish; fi- 
nally, the R3 and R5 regions merge with each other. In 
this case, the value of J2/J1 at the boundary between 
the R5 and R 6 regions decreases across J2/J1 = 1- In 
any cases, it is important to note that from our finding 
of an intermediate phase in all the three cases in Fig. Q] 
between the ferrimagnetic phase and the nonmagnetic 
phase, the phase is considered to exist irrespective of the 
strip width. 

Finally, one should note that the intermediate phase 
between the Lieb-Mattis ferrimagnetic and nonmagnetic 
states is observed in other cases. However, this phase 
is not always ferrimagnetic. One of the cases observed is 
the case of the three-leg ladder system forming a strip lat- 
tice obtained by cutting off from the spatially anisotropic 
triangular lattice in ref. IVIl in which the properties of 
the intermediate phase were unclear. Sakai's unpublished 
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FIG. 11: Magnetization curve in the R4 region. Panel (a) is 
obtained from the DMRG calculations of J2/J1 = 1 for N = 
56 (triangle) and N = 112 (square). Panel (b) is obtained by 
zooming the region near h — in panel (a). 



study suggests that the intermediate phase is nematic 38 . 
Careful examinations are required to investigate such an 
intermediate phase if it is found. 



VI. SUMMARY 

We have studied the ground-state properties of the 
5 = 1/2 Hcisenberg models on the kagome strip lattices 
depicted in Figs. [TJb) and Hfc) by the ED and DMRG 
methods. As a common phenomenon in the ground state 
of both cases, we have confirmed the existence of the non- 
Lieb-Mattis ferrimagnetism between the Lieb-Mattis fer- 
rimagnetic phase and the nonmagnetic phase. We have 
clearly found incommensurate modulations with long- 
distance periodicity in the non-Lieb-Mattis ferrimagnetic 
state. The occurrence of the non-Lieb-Mattis ferrimag- 
netism irrespective of strip width strongly suggests that 
the intermediate state found in the case of the spatially 
anisotropic kagome lattice in two dimensions is the non- 
Lieb-Mattis ferrimagnetism. 
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